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Abstract. We observed the nearby galaxy M 51 (NGC 5194) with BeppoSAX. 
The X-ray properties of the nucleus below 10 keV are almost the same as the 
ASCA results regarding the hard component and the neutral Fe Ka line, but the 
intensity is about half of the ASCA 1993 data. Beyond this, in the BeppoSAX PDS 
data, we detected a bright hard X-ray emission component which dominates above 
10 keV. The 10 - 100 keV flux and luminosity of this component are respectively 
2 X 10"^^ erg s~^ cm" and 2 x 10*^ erg s"\ These are about 10 times higher than 
the extrapolation from the soft X-ray band, and similar to the flux observed with 
Ginga, which found a bright power law component in 2 - 20 keV band. Considering 
other wavelength properties and the X-ray luminosity, together with strong neutral 
Fe K line, the hard X-ray emission most likely arises from a low luminosity active 
nucleus, which is obscured with a column density of ~ lO^^cm"^. This suggests 
that hidden low luminosity AGNs may well be present in other nearby galaxies. We 
interpret the discrepancy between Ginga and other X-ray satellites to be due to a 
large variability of absorption column density toward the line of sight over several 
years, suggesting that the Compton thick absorption material may be present on a 
spatial scale of a parsec. Apart from the nucleus, several ultra-luminous off-nuclear 
X-ray sources detected in M 51 exhibit long-term time variability, suggesting the 
state transition similar to that observed in Galactic black hole candidates. 
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1. Introduction 



M 51 is a face-on spiral galaxy, also known as the Whirlpool galaxy; it is located at 



a distance of 9.6 Mpc (Sandage & Tammann 1975 ). Optical studies of emission lines 
classified it as a LINER (Low Ionization Nuclear Emission Region) or a Seyfert 2 galaxy 



(Stauffer 1982; Filippenko & Sargent 1985; Ho et al. 1997b). Furthermore, Ho et al 



(1997c) suggested the presence of a broad Ha line. Kohno et al. (1996) found a nuclear 
molecular disk, and constrained the dynamical mass within 70 pc of the nucleus to be 
(4 - 7)xlO^M0. This indicates that M 51 hosts as massive a black hole as many bright 
AGNs. 

Observations in other wavelengths also strongly suggest the presence of an AGN 
in M 51, although previous X-ray observations of M 51 gave puzzling results. Einstein 



(Palumbo et al. |1985| ) and ROSAT (Marston et al. |1995| Ehle et al. |1995D observations 
constrained the soft X-ray luminosity of the M 51 nucleus to be Lx < 5 x lO'^^erg s^^. 
Ginga scanning observations in 1988 detected bright hard X-ray emission with Lx ~ 
1.2 X lO^'^erg/s (2 - 20 keV) with a photon index of 1.43 ± 0.08 (Makishima et al. |1990| ) 
and an intrinsic absorption of < 7 x lO^^cm^^. Such hard nuclear X-ray emission is 
often considered to be evidence of Low Luminosity AGN (LLAGN) , but it should have 
been detected with Einstein and ROSAT if the flux did not change significantly. After 
the Ginga observations, ASCA observed M 51 in the hard X-ray band (2-10 keV) 
in 1993, and did not detect such a bright hard component (Terashima et al. |l998a ). 
Instead, a faint hard X-ray continuum with a neutral Fe K line was detected, whose 
flux was an order of magnitude lower than that measured by Ginga. One possibility is 
that the difference between these observations is due to time variability with a large 
amplitude. Such a strong time variability is not often observed for low luminosity AGNs, 
and therefore is quite interesting if it is correct. The other possibility is a strong change of 
the intrinsic absorption. The key phenomenon that can help in resolving this mystery is 



the presence of a neutral Fe K line emission detected by ASCA (Terashima et al. |l998a| ), 
which indicates that during the ASCA observations, the nucleus of M 51 was heavily 
obscured by Compton thick material below 10 keV, as found in some Seyfert 2 galaxies 



(e.g. Risaliti et al. 199£ ). This idea can be confirmed by X-ray observations at energies 
above 10 keV. 

Detection of an obscured hard X-ray emission from M 51 would suggest that many 
galaxies may harbor previously unknown hidden LLAGNs emitting below 10 keV. 
Moreover, such a detection can be valuable for a study of the physical structure of 
Compton thick absorbing material around AGNs. To further investigate such a possibil- 
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ity, we performed an X-ray observation of M 51 with BeppoSAX, sensitive to hard X-ray 
emission (Boeha et al. 1997b). 



2. Observations and Data Reduction 

We observed M 51 with BeppoSAX on Jan. 18-20, 2000. The three Narrow Field 



Instruments, LEGS (Parmar et al. |l997j ), MECS2+3 (Boella et al. |l997bD , and PDS 



(Frontera et al. 1997| ) collected the data during the observation. We obtained the net 



exposures of 37 ks, 97 ks, and 45 ks, for LEGS, MEGS, and PDS, respectively. We uti- 
lized the cleaned and linearized data supplied by the BeppoSAX Science Data Center. 
The light curves, spectra, and images were reduced with the standard software, such as 
XSELEGT and FTOOLS-4.2. Throughout this paper, the data of two MEGS instruments 
were summed together. The cross-normalization constant between MEGS and PDS is set 
to be 0.85 in this paper (Fiore et al. 1999| ). 



In addition, we utilized ASGA (Tanaka et al. 1994) archival data in order to investigate 
the time variability of the M 51 nucleus. ASGA observed M 51 twice, on May 11-13, 1993 
and May 5-7, 1994. The results of ASGA 1993 observation were reported by Terashima 



et al. (1998a). Ptak et al. ( 1999] ) also performed the analysis of both 1993 and 1994 data. 



but they did not consider the contamination by the neighboring off-nuclear bright X-ray 
source. The screened data were obtained from the ASGA archival data center. In the data 
reduction, we imposed the selection criterion on the GIS data requiring that the rigidity 
is > 6 GeV c^i (Ohashi et al. |1996| ; Makishima et al. [1996D . The net GIS exposure is 



35-38 ksec for both observations, but it was only 12 ksec for the SIS data in 1994 due to 
problems of the onboard data processing. 

3. Image Analysis 

Before spectral analysis can be performed on the nuclear component, it is necessary to 
investigate X-ray images because M 51 is known to contain several non-nuclear X-ray 
sources. In Fig. 1, we show the X-ray images taken with the BeppoSAX MEGS. The 
MEGS image in the total band is overlaid on the optical image taken from SkyVIEW 
from NASA/GSFG. Three X-ray sources are clearly visible: the nucleus of M 51, a source 
in the companion galaxy NGG 5195 (source A), and a source to the east of the nucleus 
(source B). The nucleus is particularly bright in the 5.7 - 7.0 keV band, where a neutral 



Fe K line exists. Recently, Terashima et al. ( 2001 ) confirmed that this source is resolved 



as a nuclear point source with Ghandra. The source B is not the supernova SN1994I, 
since its position is different and the flux observed in ROSAT is quite low (Immler et al. 
199^ ). 



In Fig. 2, we show ASGA GIS images of M 51 taken in 1993 and 1994. The nucleus 
and the source A are always visible, in the ASGA 1993 and BeppoSAX 2000 observations. 
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The source B is not seen in the ASCA GIS images, indicating the time variabihty of this 
source. On the other hand, the north-east X-ray source (source C) that was visible in the 
ASCA 1994 observation is not seen in the ASCA 1993 and BeppoSAX image. 

Sources B and C are thought to correspond to the ROSAT PSPC point sources R3 and 



Rl in Marston et al. (1995), respectively. All the detected sources have X-ray luminosities 



in the range of ^ 10'^^ erg s ^. 

4. Spectral Analysis 

4.1. BeppoSAX Spectra of the Nucleus 

In order to investigate the spectrum of the nucleus, we extracted data from the imaging 
instruments onboard BeppoSAX corresponding to a radius of 2' so that neighboring 
X-ray sources do not contaminate the spectra of the nuclear source. The background 
spectra were prepared using data from the same detector region, utilizing the standard 
background data set (Nov. 1998 for MECS, Dec. 1999 for LECS), and we subtracted 
them from the on-source spectra. In the spectral fitting, we used the standard response 
matrices released in September 1997. 

In Fig. 3, we show the X-ray spectra of LECS and MECS. The data cannot be fit by a 
simple power law; they reveal both Fe-L and Fe K line features, but the Raymond-Smith 



thermal model (Raymond & Smith 1977 ) with the Galactic absorption of 1.5 x 10^°cm ^ 



(Stark et al. |l992D does not fit the data either: the reduced chi-square is ~ 2 for both 



models. These spectral features are similar to the ASCA 1993 spectra. Thus, we fit the 
spectra with a model including Raymond-Smith thermal plasma plus a power law, and 
including an additional narrow line. The metal abundances are fixed to 0.2 x Solar with 



the Solar metal abundance ratios (Anders & Grevesse 1989). Such a model was previously 



applied to the ASCA 1993 spectra (Terashima et al. 1998a). The fit is acceptable as shown 
in Fig. 3 and Table 1. The power law photon index of 1.59 ± 0.23 is not well constrained, 
but it is consistent with the ASCA results. The Fe K line center energy is 6.47 ± 0.14 
keV, indicating that it is due to fiuorescence from the neutral iron, in agreement with the 
ASCA results. The equivalent width of the Fe K line is 0.9^q'5 keV, which is within the 
error range of the 1993 ASCA results. The temperature of the Raymond-Smith plasma is 
0.51 ±0.16 keV. The flux and luminosity of the power law component are 3.4]';q 3 x lO""*^^ 
erg s^^ cm^^ and 3. Giro's x 10'^^ erg s~^, respectively in 2 - 10 keV band, about 30% of 
1993 ASCA results. The strength of Fe K hue is also smaller than that of 1993 ASCA 
value of (4 - 14) x lO"*" c s"^ cm"^. 

We conclude that the spectral properties below 10 kcV measured in the BeppoSAX 
data are almost the same as the ASCA results except for the flux. To extend the bandpass 
beyond that of the imaging instruments (10 keV), we plot a PDS spectrum together with 
the LECS and MECS spectra in Fig. 4a. The solid line represents the best-fit model 
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Table 1. Results of joint spectral fittings of M 51 nucleus with various models. 
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LEGS + MEGS 

LEGS + MEGS 

0.51 ±0.16 4.4+^;^ 


1.59 ± 0.24 


j^-^ ^ + 0-29 
-^-^■-^-0.22 

7 1 + 2.6 
' ■'--2.1 


ABS''*PO'' 

ABS=*(RS+PO° + GA) 
6.47 ±0.14 3.4 ±1.5 




2.19 
1.18 


LEGS + MEGS + PDS 


1.59 ± 0.23 


' -"-2.0 


ABS''*(RS+PO'+GA+ABS''*PO'') 

6.47 ±0.14 3.4 ±1.5 lot^^ 3.3^1° 




1.14 



LEGS + MEGS + PDS ABS^*(RS+PO^ + GA+ABS''*PO'') 

0.52+^-^^ 4.40+^-| 1.60 ±0.23 '^■'^tl'l 6.47 ± 0.14 3.4 ± 1.5 ^-^tl'l 1.9 fix 4.9+^;^ 1.21 

ABS^ and ABS'^: Photoelectric absorption of the Galactic and intrinsic material, RS: Raymond- Smith thermal model, PO'^ 
and PO'^: Power law model, GA: Gaussian line. 
g: Parameters of a gaussian model. 

The unit of normalizations (norm) of power law and gaussian is c s~^ cm" keV"^ at 1 keV and c s"^ cm , respectively. 



obtained by the LEGS and MEGS spectral analysis. It can be seen that the PDS flux is 
higher by a factor of ^ 10 compared with the extrapolation of the best-fit model of LEGS 
and MEGS, and significant photon flux is detected up to 50 keV. This spectral feature is 
quite similar to that seen in Gompton-thick Seyfert 2 galaxies such as NGG 4945 (Done 



et al. [1996D , Gircinus Galaxy (Matt et al. |1999|), and Mrk 3 (Gappi et al. |1999|) 

We add an absorbed power law component to represent this hard X-ray emission, 
and fit the LEGS -I- MEGS -I- PDS spectra simultaneously. The best-fit results are shown 
in Table 1 and Fig. 4b. The wide-band X-ray spectra require a hard X-ray component 
strongly absorbed by optically thick material with a column density ~ 10^^ cm~^. The 
best-fit parameters other than this newly included component do not change significantly. 
Since the photon index and absorption column density for the hard X-ray component are 
not well constrained separately (as they are correlated in spectral fits), to estimate the 
absorbing column we fixed the photon index at 1.9, typical for Seyfert 1 galaxies without a 
reflection component. Then, the absorption column density becomes 5.6l^^'g x 10^'* cm"^, 
and the X-ray flux and luminosity are 1.9 x 10^^^ erg s~^ cm"^ and 2.0 x 10''"'^ erg s"^, 
respectively in 10 - 100 keV band. The absorption - corrected X-ray flux is 2.2 x 10^^^ 
erg s~^ cm~^ in 2 - 20 keV band, similar to the Ginga results, and 1.5 x 10~^^ erg 
s~^ cm"^ in 2 - 10 keV which should be easily detected with ASGA if no absorption 
was present. When we fit only the PDS spectrum with an absorbed power law model 
by fixing the power law photon index to be 1.9, we obtained only the upper limit of 
A'h < 5.4 x lO^^'cm"^. This is most likely because the model adopted by us, with a 
luminous but heavily absorbed hard component plus a soft component, has an apparent 
absorption feature located in the energy band between MEGS and PDS. To investigate a 
possible presence of a Gompton reflection component, we include it by fitting the data to 
the spectral model "pexrav" as implemented in XSPEG; with this, the reflection fraction 
becomes <18%. 
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4.2. ASCA spectra of the Nucleus 



Since the result of the 1993 ASCA observation are reported in Terashima et al. ( 1998a ), 
we here show the result of the 1994 observation. Since the exposure time of the SIS 
is short, we created only the GIS spectrum. The source C is located at ^ 2' from the 
nucleus. Thus we created spectra using counts detected only within 1.5' of the nucleus, 
where the source C may contaminate the nuclear spectrum by ~20%. The spectral model 
used for the fits of the ASCA data is the same as that for LEGS and MEGS spectra: 
Raymond-Smith plasma, plus power law, and one emission line. The normalizations of 
thermal and power law components are left free, but the other model parameters are 
fixed to the values measured in the 1993 observation, due to the poor quality of 1994 
data. 

The spectrum and best-fit model are shown in Fig. 5, which tells us that the overall 
property of the spectrum apart from the normalization is almost the same as that of 
ASCA 1993 data and BeppoSAX data. The X-ray flux of the power law component is 
- 3 X IQ-^^ erg s^^ cni"^ jn 2 - 10 keV band, about 30% of that measured in the 1993 
data and almost the same as BeppoSAX data. The details of the Fe K line are ambiguous 
due to the small extraction area. . In order to constrain the Fe K line strength better, we 
extended the integration radius to 3.75'. Taking into account that about 45% of photons 
in this spectrum come from the source C, we obtain an upper limit on the Fe K line 
strength of 2.7 x 10^^ photons s^^ cm^^, significantly smaller than that measured in the 
1993 ASCA data. 



4.3. Time Variability of the Nucleus 

We investigate the short time variability properties of the hard component of the nucleus. 
We make light curves in the 3-10 keV band (MEGS) and 10 - 100 keV band (PDS). 
The MEGS light curve shows a counting rate from the region within 2' of the nucleus 
without background subtraction, since the count rate of the background is only 25%. The 
PDS light curve is created after background subtraction. We do not find any significant 
short-term time variability in either light curve. 

In Fig. 6, we show the long-term time variability of the X-ray flux of the power law 
component in the 2-10 keV band and absorption-corrected 2-20 keV band. The X-ray 
flux in 2 - 10 keV band represents the scattered component of the nuclear emission, and 
in 1994, it drops to 30% of the 1993 value. Its flux in the 2000 BeppoSAX data is almost 
the same as that of 1994 ASCA observation. The absorption-corrected flux in 2 - 20 keV 
of Ginga and BeppoSAX PDS corresponds to the direct nuclear emission. Nominally the 
2000 BeppoSAX flux is twice the 1988 Ginga value, but the errors are large, and these 
fluxes are consistent with the values being the same. The Fe line equivalent width is 
constant within the error range as described in §4.2. 
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4.4. Off-nuclear X-ray Sources 



M 51 includes several off-nuclear point- like sources. The typical X-ray luminosity of the 
brightest X-ray sources is lO-^^"^*^ erg s~^, much higher than the Eddington limit for 
accreting neutron stars. In the ASCA and BeppoSAX observations, three off-nuclear 
sources were visible (Figs. 1 and 2). The source A is associated with NGC 5195, and its 
flux appears constant. The other two sources (B and C) are located in the spiral arms. 
They exhibit strong time variability, and are visible when they become as bright as the 
nucleus. 

The source B shows a large time variability in 2000 as shown in Fig. 7, where the 
count rate in the second half of the observation decreases to ~50% of that measured in 
the first half. This is the third example of short-term time variability of the off-nuclear 
ultra-luminous X-ray sources, after those in IC342 (Okada et al. |l998| ; Kubota et al. |200l|) 



and M82 (Matsumoto et al. 1999 ). The spectra of the sources A and B are determined 
from the counts collected within 2' of their position. For the source A, we prepared the 
spectra using the ASCA 1993 data and the BeppoSAX data. We do not use the ASCA 
1994 data because the CIS instrument cannot resolve the source A clearly (see Fig. 2). 
The source B is only visible in the BeppoSAX data, and its count rate varies significantly. 
We thus analyzed two MFCS spectra; one is integrated in the bright phase (0 - 110,000 
sec in Fig. 7), and the other is in the faint phase (110,000 - 190,000 sec in Fig. 7). 

The X-ray spectra of these ultra-luminous off-nuclear X-ray sources are known to be 



generally described by multi-color disk black body (MCD: Mitsuda et al. 1984; Makishima 
et al. |198e| ; Mizuno et al. |1999| ; Makishima et al. |2000| ; Kotoku et al. pOOl] ), and thus we 
try to fit the spectra by a power law model and MCD model. The spectrum of source A 
is almost the same between the 1993 ASCA and BeppoSAX observations, and thus we 
fitted both data sets simultaneously. We performed spectral fittings with the energy band 
above 2 keV, to avoid the uncertainty of the diffuse thermal emission due to the different 
angular resolution among instruments. Instead, we include the Raymond-Smith thermal 
model to represent the diffuse emission. Its temperature cannot be well determined, and 
thus we fixed the temperature of Raymond-Smith model to 0.39 and 0.35 keV for source 
A and B, respectively, based on the spectral fitting of the ASCA SIS data around each 
source. Since we are interested in the hard component rather than soft thermal emission, 
we performed the spectral fitting above 2 keV. The results are summarized in Table 2 
and Fig. 8. 

The limited photon statistics does not allow us to distinguish between these two 
models for both sources. The inner disk temperature of MCD model for the source A, 
1.5 - 2.2 keV, is typical for ultra-luminous off-nuclear X-ray sources in other galaxies. 
The luminosity of the MCD component in the 2-10 keV band is 3.2 x lO'^^ erg s~^ 
for source A. For the source B, the spectrum in the high flux period seems to be harder 
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Table 2. Results of spectral fittings of source A and B with various models above 2 keV. 
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r: parameters of Raymond-Smith model 
p: parameters of Power law model 
d: parameters of MGD model 

The unit of norm*' is 10"^c s"^ cm"^ keV"^ at 1 keV. 

The norm'' represents ((iiin/lkm)/(_D/lOkpc))^cos0 where Rin and 6 are an inner radius 
and an inclination angle of the disk, respectively, and D is a distance to the source. 



than that in the low flux period; the power law photon index or the temperature of the 
MCD model is significantly different. We expect a constant radius of the MCD model 
despite different temperatures, but that cannot be tested due to the poor data quality. 
The luminosities of the MCD component in the 2-10 keV band of the source B are 
3.8 X lO'^^ erg s~^ and 1.2 x 10"^^ erg s~^ for the bright and faint state, respectively. 

In order to study the long-term time variability of these sources, we utilized the 
archival ROSAT HRI data from HEASARC (NASA/GSFC High Energy Astrophysics 



Science Archive Research Center) and the paper of Ehle et al. ( 1995 ). We cannot resolve 
the source A clearly in the HRI images, and thus the source A may be an assembly of 
discrete sources in NGC 5195. The constant flux is consistent with this scenario. The flux 
histories for the sources B and C are summarized in Table 3. These are typically in the 
range of (0 - 3) x 10~^^ erg cm^^ s^^, at a level comparable to the ASCA and BeppoSAX 
measurements. The source B varies by a factor of 2, while the variability of the source C 
is quite strong, by a factor of >5. 
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Table 3. ROSAT HRI Flux of the Source B and C 



Time 



1991.12.07 



1992.01.08 
1992.05.22 - 25 
1992.06.05 



1994.05.22 - 23 1994.06.18 - 24 1997.12.26 - 30 



exposure 16097s 
B 1.5 
C 2.2 



10411s 
1.2 
<1.0 



9392s 
1.8 
2.4 



36343s 
2.1 
1.7 



8104s 
0.9 
<0.4 



The unit is 10 ^'^ erg s ^ cm ^ . The conversion from counts to source flux density is 4.08 x 

-11 



10" 



erg cm ^ count ^ , from Ehle et al. (1995). The data in 1992 are quoted from Ehle et 



al. (1995). 



5. Discussion 

5.1. Bright Hard X-ray Emission 

In our observations of M 51 with BeppoSAX, we detected in the PDS data bright hard 
X-ray emission above 10 keV, which we interpret as associated with an active nucleus 
present in the galaxy. This inference is supported by the fact that in the MECS data we 
detect a neutral Fe K line emission from the central point-like source, and such emission 
is consistent with the ASCA observations. The flux above 10 keV is about 10 times 
higher than the extrapolation of the nuclear flux detected by the MECS below 10 keV. 
Considering these spectral features, we infer that M 51 harbors a heavily obscured Seyfert 
2 type nucleus, with the column density 10^'* cm~^ 

The X-ray luminosity of this hard component is 2.0 x lO''^ erg s~^ in the 10 - 100 
keV band, and thus M 51 is probably the first low-luminosity Conipton thick Seyfert 
2 AGN. The existence of AGN with X-ray luminosity of ~ lO"'^ erg s~^ is reasonable, 
according to the prediction from the correlation of Ha and X-ray luminosity in Seyfert 1 



galaxies (Ward et al. |1988), and the correlation of [Oni] and X-ray luminosity of Seyfert 



galaxies (Mulchaey et al. 1994), based on the M 51 value of Lqiii ^ -^Hq ~ 10'^^ erg s~^. 
Moreover, the presence of an obscured nucleus in M 51 is suggested by a dense molecular 
disk with > 3 x lO^'^cm"^ detected around the nucleus by HCN and CO observations 



(Kohno et al. 1996) 



Nevertheless, we cannot completely rule out the possibility that the PDS spectrum 
is contaminated by bright hard sources, because the absorption is ambiguous if one 
considers the PDS spectrum alone. However, there are several observational results which 
support the above conclusion. The Ginga observation - using the scanning mode, and 
thus allowing good positional accuracy - indicated that the bright hard X-ray emission 



indeed comes from the M 51 nucleus (Makishima et al. 1990 ). Since the hard X-ray flux 



is similar between the Ginga LAC and the BeppoSAX PDS, the hard emission detected 
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from the PDS is also likely to come from M 51. In addition, there are no other bright 
X-ray point sources with a flux > 10^^^ erg s~^ cm~^ (0.1 - 2 keV) within 2 degrees of 
the M 51 nucleus in the ROSAT PSPC image. Consequently, we conclude that the hard 
X-ray emission detected in the PDS is associated with the nucleus of M 51. Then, M 51 
is found to be one of the brightest low luminosity AGNs in the sky above the 10 keV 



band, together with M 81 and NGC 4258 (Pellegrini et al. |2000| ; Fiore et al. |200l|). It is 
interesting that these two objects are not Compton thick, unlike M51. 

There still remains a difference between the Ginga and other satellite observations; 
the X-ray flux below 10 keV obtained with Ginga was higher than that measured by other 
instruments by a factor of more than 20, and the X-ray spectrum obtained with Ginga 
exhibited no absorption features. One possibility is the change of the column density of 
the absorbing material in the line of sight; the absorption happened to be weak in the 
Ginga observation. Such a case is in fact observed for the famous Seyfert galaxy NGC 



4051 (Uttley et al. |l999|) , and for NGC 1365 (Risaliti et al. |2000|) with BeppoSAX, and 
thus such a behavior is not unique to M 51. Another scenario would have a large change 
in the soft (< 10 keV) flux, by a factor of more than 20. However, the X-ray flux above 10 
keV band is almost the same between the Ginga and BeppoSAX observation epochs, and 
in addition agrees with the prediction from Lqiii a-nd Lhq, suggesting that the former 
explanation is more likely. 



5.2. A Neutral Fe K line and Reflected Continuum 

BeppoSAX detected neutral Fe K lines from the M 51 nucleus, confirming the previous 



ASCA detection (Terashima et al. 1998a ). This line is likely to arise via fluorescence 
from neutral material. The equivalent width of the Fe K line, 0.9^q'5 keV is consistent 



with the ASCA 1993 result (Terashima et al. |l998aj). The equivalent width of the Fe 
K line can be explained as arising via reflection by the reprocessing material covering 
27r around the nucleus with Ix Solar metal abundances. Other well-known examples 



of objects indicating such geometry are NGC 6552 (Fukazawa et al. 1994) and Circinus 



Galaxy (Matt et al. 1996) 



A natural association of the Fe K line would be with the dense, neutral material 
present around the nucleus; such material is most likely responsible for the Compton 
reflection component, seen in many Seyfert 1 galaxies. We investigated how strong of 
such reflection component is in the BeppoSAX data, but obtained a small upper limit 
corresponding to a solid angle subtended by the reflector of 18% of 27r. If the material 
responsible for the Fe K line and the Compton reflection were to be the same, this value 
seems to be in contradiction with the Fe K equivalent width. A more likely scenario is 
that the Fe K line is produced in a Compton-thick region in transmission. In fact, the 
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equivalent width of the Fe K line can be reproduced by such a situation with a column 
density of ^ lO^'^^cm'^, covering factor of 2tt, and a Solar Fe abundance (Makishima 198£). 



5.3. Time variability 

The bright hard component that was detected with Ginga and BeppoSAX, presumably 
due to direct nuclear emission, has almost constant flux during the BeppoSAX observa- 
tion. Such a trend of a weak time variability is similar to that seen in other LLAGNs 
(Ishisaki et al. |l996t lyomoto et al. |1996| , ^00l| lyomoto |l999t Reynolds et al. |200C| ), and 



is not consistent with the trend seen in Seyfert galaxies (Nandra et al. 1997; Hayashida 
et al. |1998| ). We note that a scenario involving an advection dominated accretion flow 
(ADAF) is consistent with such a behavior, but any more detailed discussion of the ap- 
phcability of ADAFs to LLAGN is beyond the scope of this paper (but see, e.g., Lasota 



et al. [1996D . 

On the other hand, we detected time variability of two quantities: the weak hard X-ray 
emission below 10 keV, and the absorption column density for direct nuclear emission. 
These are indicative of a change of conditions of the material surrounding the nucleus. 
The former, which is thought to be a scattered or reflected nuclear component, exhibits 
a long time variability, by a factor of 2 - 3 on a time scale of years, through ASCA and 
BeppoSAX observations (Fig. 7). The strength of the Fe K line also varied simultane- 
ously with the continuum by almost the same factor. Regarding the latter, we measure 
the intrinsic absorption column density changes by a factor of > 100 within several 
years between the Ginga and ASCA observations. These indicate that the Compton- 
thick absorption material has a spatial scale (meaning the size and/or a distance from 
the nucleus) of a parsec. Such a parsec - scale torus material is also suggested by other 
observations (Greenhill et al. 1996| ; Gallimore et al. 1997, Siebenmorgen et al. 1997; Matt 



200C ) . Further variability studies of direct nuclear emission on shorter time scales at en- 



ergies above 10 kcV, such as with INTEGRAL and Astro-EII HXD, are necessary to 
characterize the dimensions of this component in more detail. 



5.4. Obscured Low Luminosity AGNs 

ASCA survey showed that LINER galaxies often exhibit hard X-ray emission with X-ray 
luminosity of (2 - 5) x lO-^Oerg/s (2-10 keV) and power law spectra with photon indices 



of -1.7 (Makishima et al. |1994[ Ishisaki et al. |l996t lyomoto et al. |1996| ; |1997| ; |1998| ; 
Terashima et al. 1998a; 1998b; pOOO ). However, there are LINERs which do not show 
a significant excess of hard X-ray emission below 10 keV. Interestingly, ASCA results 
show that some LLAGNs also exhibit an intrinsic absorption of — 10^'^cm^^ which is 
comparable to Compton-thin Seyfert 2 galaxies. Furthermore, ASCA found two LLAGNs, 



M 51 and NGC 1365 (Terashima et al. |1998a| ; lyomoto et al. |1997| ), which exhibit a 
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flat hard X-ray component with a strong neutral Fe K line, and thus these are good 
candidates for the low luminosity Compton-thick Seyfert 2 galaxies. In this context, our 
BeppoSAX observations confirm M 51 as a Compton-thick LLAGN. These results suggest 
that the distribution of the absorption column density of LLAGNs may be similar to that 
of Seyfert galaxies, and thus a significant fraction of LINERs may be also completely 
obscured below 10 keV, as is the case for Seyfert galaxies. LLAGNs are thought to be 



quite numerous because LINERs are quite numerous (Keel 1983; Ho et al. 1997a). With 
this, future X-ray observations of LINER galaxies may reveal many LLAGNs in galaxies 
that are considered otherwise "normal." Understanding the relationship between LINERs 
and LLAGN is likely to be an important step in understanding the link between evolution 
of normal and active galaxies. 



5.5. Ultra-Luminous Off-nuclear X-ray Sources 



M 51 contains several ultra- luminous off- nuclear X-ray sources found with Einstein 
and ROSAT (Palumbo et al. |l985[ Marston et al. [l995| ; Ehle et al. |l995| ). ASCA and 
BeppoSAX observations resolved three of them, and all three exhibit time variability of a 
factor of ^ 3. The X-ray luminosity measured in those sources - 10-^^"*° erg s~^ - cannot 
be explained by neutron star binaries, as the Eddington limit for those is 2 x 10'^^ erg 
s~^. The X-ray spectra for two of them can be fitted with power law or disk blackbody 
models. This spectral shape is similar to that seen in other ultra-luminous off-nuclear 
X-ray sources in other galaxies (Mizuno et al. |1999| Makishima et al. |2000| ; Mizuno [2000| ) . 
These bright X-ray sources in nearby galaxies are postulated to be black hole binaries, 
where spectra (and in particular, the disk black body temperature and normalization) 
imply Kerr black holes with masses -100 Mq (Mizuno et al. |1999| Makishima et al. |200(]| ; 



Mizuno 2000). We suggest that the bright, off- nuclear X-ray sources in M 51 detected 



with ASCA and BeppoSAX are objects of the same type. The strong variability of these 
sources is similar to that often observed in the Galactic black hole candidates. Since their 
X-ray spectra were obtained only in their bright states, it would be quite interesting to 
know what are the X-ray spectra in their faint states, and in particular, whether they 
exhibit hard power law spectra in their faint states, similar to the Galactic black hole 
candidates. So far, two ultra-luminous off nuclear X-ray sources in IC342 exhibited such 



a spectral transition (Kubota et al. 2001 ), suggesting that these are black-hole binaries. 
Since it contains relatively many ultra- luminous off- nuclear X-ray sources, M 51 is an 
interesting target for repeated observations by Chandra and XMM-Newton to study such 
sources. 

We are grateful to Prof. Greg Madejski for correction of the English in the manuscript. 
The authors also thank the anonymous referee for careful reading and helpful comments, 
and thank as well operation and calibration teams of BeppoSAX and ASCA. 
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Fig. 1. BeppoSAX MECS images of M 51. (a) MECS image (contour) in the total band 
overlaid on the optical image (gray scale). (b)MECS image in the energy band of 5-7 
keV. The center position is (13h 30m 00s, 47d 13m 00s). The axis unit is arcmin. 
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Fig. 2. ASCA GIS images of M 51 in 1993 (left) and 1994 (right). The axis unit is arcmin. 
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Fig. 3. LEGS and MEGS spectra of the nueleus. The sohd hne represents the best-fit 
model (Raymond-Smith model -|- power law model -|- Gaussian). 




Fig. 4. LEGS, MEGS, and PDS spectra of the nucleus. The solid line represents the best- 
fit model (the same as Fig. 4) for the LEGS and MEGS spectra (left), and the best-fit 
model (Raymond-Smith model -|- power law model -|- Gaussian -|- absorbed power law) 
for all the instruments (right). 
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Fig. 5. ASCA GIS spectrum of the nucleus in 1994. The solid line represents the best-fit 
model (see text). 
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Fig. 6. Light curve of the M 51 nucleus measured by various satellites. The open circles 
represent the flux of the weak continuum in 2 10 kcV band, which arc thought to be 
a scattered or reflected component. The open triangles are the flux of the bright nuclear 
direct continuum in 2 - 20 keV band. 
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Fig. 7. MECS light curve of the ofF-nuclear source B. 
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Fig. 8. X-ray spectra of the source A and B. (a)ASCA GIS/SIS (1993) and BeppoSAX 
LECS/MECS (2000) spectra of the source A. (b)MECS spectra of the source B. The left 
and right are obtained in the bright and faint phase, respectively. 



